Tri-partite Entanglement Witnesses and Sudden Death 
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I explore entanglement dynamics in a three qubit system comparing the ability of entangle- 
ment witnesses to detect tri-partite entanglement to the phenomenon of entanglement sudden death 
(ESD). Using a system subject to dephasing I invoke entanglement witnesses to detect tri-partite 
GHZ and W-type entanglement and compare the evolution of their detection capabilites with the 
evolution of the negativity, bi-partite concurrence, and tri-partite negativity. Interestingly, I find a 
state in which there is no concurrence or tri-partite negativity but there is entanglement. Finally, I 
utilize a three qubit quantum error correction (QEC) code to address how ESD affects the abilities 
of quantum error correction. 

PACS numbers: 03.67.Mn, 03.67.Bg, 03.67.Pp 



I. INTRODUCTION 

Entanglement is a uniquely quantum mechanical phe- 
nomenon in which quantum systems exhibit correlations 
not possible for classical systems. As such, entanglement 
is a vital resource for many aspect of quantum informa- 
tion processing including quantum computation, quan- 
tum metrology, and quantum communication [l[ . But de- 
spite its fundamental and practical importance and much 
work in the subject, there are many aspects of entangle- 
ment, especially multi-partite entanglement, that are in 
need of further study Q • 

A major challenge facing experimental implementa- 
tions of quantum computation, sensing, and communica- 
tion is decohcrence, unwanted interactions between the 
system and environment. Decoherence may be especially 
detrimental to highly non-classical, and hence most po- 
tentially useful, entangled states Q. A manifestation of 
this is entanglement suddent death (ESD) in which en- 
tanglement is completely lost in a finite time 0, Q de- 
spite the fact that the coherence loss of the system is 
asymptotic. This aspect of entanglement has been well 
explored in the case of bi-partite systems and there are 
a number of studies looking at ESD in multi-partite sys- 
tems d, 0i H ■ I n addition, there have been several ini- 
tial experimental studies of this phenomenon [To| . How- 
ever, even when analyzing a multi-partite system, pre- 
vious works demonstrate ESD only for bi-partite entan- 
glement, either via concurrence or negativity rather than 
using measures for purely multi-partite entanglement. It 
is important to note that the characterization and quan- 
tification of true multi-partite entanglement is still very 
much an unsettled area for pure states and even more so 
for mixed states. 

In this paper I explore the loss of detectable entangle- 
ment in a three qubit system by invoking entanglement 
witnesses, observables that can detect the presence of en- 
tanglement. Entanglement may be present in a system 
but still not be practically useful 0]. For entanglement 
to be useful its presence should be efficiently detectable 
experimentally. Multi-partite entanglement can be de- 
tected inefficiently via quantum state tomography or a 



violation of Bell inequalities. It can be detected efficiently 
by utilizing properly constructed entanglement witnesses 
[14| . I compare the entanglement detection abilites of 
tri-partite entanglement witnesses to ESD of tri- and bi- 
partite entanglement in a given system. In this explo- 
ration I find a state which has no concurrence and no 
tri-partite entanglement as measured by the tri-partite 
negativity but is entangled as measured by the negativ- 
ity. I then apply these results to a three qubit quantum 
error correction (QEC) code and explore how ESD affects 
the working of this code. 

Three qubit pure states can assume a 'standard' form 
with respect to local unitary operations [l5| . 

|V> = A o |000) + Aie l9 |100) + A 2 |101) + A 3 |H0) + A 4 |lll), 

(1) 

where Xi > 0, ^ i Af = 1 and < < n. These states can 
be separated into four broad categories: separable (in all 
three qubits), biseparable, and there exist two types of 
locally inequivalcnt tri-partite entanglement (GHZ and 
W-type) [13]. Similar classification schemes exist for 
mixed states. Reference defines four classes of three 
qubit mixed states each of which includes the preceeding 
classes as special cases. They are separable (S) states, bi- 
separable (B) states, W states, and GHZ states, which 
encompasses the complete set of three qubit states. Note 
that additional subtlety exists in characterizing the en- 
tanglement within each of these classes [l2j ■ 

To determine in which class a given state belongs one 
can use entanglement witnesses, observables which give 
a positive or zero expectation value for all states of a 
given class and negative expectation values for at least 
one state in a higher (i.e. more inclusive) class. Specif- 
ically, I will make use of entanglement wintcsscs to 
identify whether a state is in the GHZ\W class (i. e. a 
state in the GHZ class but not in the W class), in which 
case it certainly has GHZ type tri-partite entanglement, 
the W\B class, in which case the state certainly has true 
tri-partite entanglement either of the GHZ-type or W- 
type, or the B class in which case it is not certain that 
the state has any tri-partite entanglement. While the 
witnesses we explore may not be of the sort that can be 
implemented efficiently for experimentally determining 
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the presence of tri-partite entanglement, they are among 
the most sensitive, or finest, known witnesses. Thus, if 
these witnesses do not detect the presence of entangle- 
ment neither will any of the effeciently implcmentable 
witnesses. 

I will compare the detection ability of the witnesses to 
the evolution of the concurrence fl6} , Cjk , for measuring 
the bi-partite entanglement between qubits j and k after 
partial trace over one qubit. The concurrence of a two 
qubit state with density matrix p^ is defined as the max- 
imum of zero and A, where A = ^/\[ — \f\^ — VA3 — VXj" 
and the \ are the eigenvalues of Pjk{o- y <g>a y ; )p* k (a y ®a y : ) 
in decreasing order, where a y is the y Pauli matrix of 
qubit i. Other entanglement measures that I will look at 
are the negativity, N, for which I will use the sum of the 
absolute values of the negative eigenvalues of the partial 
transpose of the density matrix [17| with respect to one 
qubit, and the tri-partitc negativity, A^, a tri-partitc en- 
tanglement measure for mixed states which is simply the 
third root of the product of the negativities with respect 
to each of three qubits [ljj • If the negativity is the same 
when taking the partial transpose with respect to any of 
the three qubits, N = N3. 

We look at a three qubit system, with no interaction 
between the qubits, placed in a dephasing environment 
fully described by the Kraus operators 

where the dephasing parameter p can also be written in a 
time-dependent fashion p = 1 — exp(— nt). When all three 
qubits undergo dephasing we have eight Kraus operators 
each of the form Ai = (Ki®Kj®Kk) where I = 1,2, ...,8 
and i,j, k = 1, 2. 

The next three sections are dedicated to exploring the 
entanglement dynamics and detectability of three types 
of three qubit states: GHZ-type states in Section [III W- 
type states in Section IIIII and a state with GHZ-type 
tri-partitc entanglement and bi-partite entanglement in 
Section IIV1 In Section |V] we apply our results to the 
three qubit phase flip quantum error correction code and 
study the affect of ESD on the code. Section IVT1 contains 
some further observations and conclusions. 



II. GHZ-TYPE STATES 

In this section I explore the affects of dephasing on the 
following three qubit GHZ-type states, 

P=^-^-^ + 1Pg (3) 

where p G = (|000) + |111))((000| + (111|). These states 
have no bi-partite entanglement and non-zero negativity 
and tri-partite negativity equal to — min [|(1 — 5q), 0] . 
The state p is certainly in the GHZ\W class for q > | 
since it can be detected via the witness [ll[ Wg = f 1 — 



Pq. p is in at least the class W\B for q > | since it can 
be detected via the witness [ll[ Wwi = |l — Pg- 

The effect of dephasing all of the qubits of p is sim- 
ply to reduce the magnitude of the off diagonal elements 
I — ► |e _3/tt / 2 . For q = 1 the entanglement witness Ww2 
gives — ie -3Kt / 2 which decays exponentially with time. 
In other words, the state can always be detected by a 
W-state witness and does not exhibit tri-partite ESD 0| • 
However, the dephased state cannot always be detected 
as belonging to the GHZ\W class. This can be seen from 
the expectation value Tr(WGp) which is negative (thus 
detected) only as long as nt < |ln(2). Thus, we have 
a 'sudden death' type phenomenon with respect to the 
GHZ entanglement witness. Experimentally we cannot 
be sure the state exhibits GHZ-type entanglement. 

The inability to detect any tri-partite entanglement 
in the state p can occur for q < 1. In such a case we 
find Tr(Ww2p) = |(3 - 3q- 4 qe - 3Kt / 2 ) and Tr(W G p) = 
g(5— 3q— 4qe -3fct / 2 ). The behavior, as a function of time, 
for both entanglement witnesses is similar, as shown in 
Fig. [U and clearly demonstrates the inability of the wit- 
nesses to detect any entanglement at sufficiently long but 
finite time. 

Comparing this to the entanglement measures, we 
note that there is only one eigenvalue of the partially 
transposed density matrix that can be negative and so 
N = N 3 = - min [|(1 - q - Aqe~ ZKt / 2 ), 0] . For q = 1, 
i.e. initially pure states, this is equal to — Tr(Ww2/o) 
and thus the witness always detects the tri-partite en- 
tanglement. For q < 1 the initial state is mixed and 
— Tr(Ww2P) is no longer equal to N. Instead, the entan- 
glement witness goes to zero at a rate approaching three 
times that of A" demonstrating that there is remaining 
tri-partite entanglement in the system beyond the de- 
tection of the entanglement witness. Nevertheless, even 
N decays to zero in finite time exhibiting the ESD phe- 
nomenon, Fig. HJ As expected, C between any two qubits 
of p after partial trace over the third qubit is always zero 
thus the system never exhibits bi-partite entanglement. 

III. W-TYPE STATES 

Unlike GHZ states the W state \W) = -75GOOI) + 
|010) + 1 100) ) , retains a high degree of (bi-partite) entan- 
glement upon partial trace over one qubit. Thus, W-type 
states allow some comparison between tri-partite and bi- 
partite entanglement evolution. We start with the state 

p=^-± + qpw (4) 

where pw = \W)(W\. States of this sort in the W\B 
class may be detected by the entanglement witness [ll| 
WVi = |l — pw For q = 1 the expectation value of 
the entanglement witness is — ^ and the concurrence of 
any two qubits is |. The entanglement witness detects 
the state, for all q > — .619. C, after partial trace 
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FIG. 1: (Color online) Demonstration of tri-partite entan- 
glement sudden death in GHZ states. Top-left: Expectation 
value of the entanglement witness — Tr(WGp) as a function of 
the initial state paramaterized by q, and the dephasing con- 
stant Kt. At sufficiently low values of q or high values of Kt 
the state can no longer be detected by a GHZ-type entangle- 
ment witness. Top-right: — Tr(Ww2p) as a function of the 
same parameters. This entanglement witness shows when the 
state can no longer be detected as one with tri-partite entan- 
glement. Bottom: The negativity with respect to any qubit 
partition (and thus in this case the negativity is equivalent 
to the tri-partite negativity) . When this is zero the state has 
no distillable entanglement. For q — 1 (initial pure state) the 
negativity is equal to — Tr(Ww2p) meaning that the VVw2 
entanglement witness always detects the presence of the tri- 
partite entanglement. When q < 1 the state is mixed and the 
witness does not always detect the tri-partite entanglement 
which is known to be present via the tri-partite negativity 
measure. 




FIG. 2: (Color online) Demonstration of bi- and tri-partite 
entanglement sudden death in W-states. Each plot shows the 
expectation value of the entanglement witness — Tr(Wwpw) 
(solid line), the bi-partite concurrence (medium dashed line), 
and the negativity (large dashed line) as a function of the 
dephasing constant Kt. Top left: q = 1, the initial state is 
pure. In this case — Tr(Wwpw) goes to zero in finite time, but 
the other entanglement measures do not. Top right: q = .95, 
for any q < 1 the concurrence undergoes ESD before the 
negativity (which is equal to the tri-partite negativity) goes 
to zero. Bottom: q — .75. 

partite entanglement witness going to zero first followed 
by C and finally the N = N3. In other words, after a 
certain time there is no longer bi-partite entanglement in 
the system despite the remaining W-typc tri-partite en- 
tanglement. The above expressions also show that each 
entanglement measure decreases at the exponential rate 
— Kt and there is no difference in the exponential rate 
of decay between the bi-partite and tri-partite entangle- 
ment. The difference in entanglement measure evolution 
comes from different coefficients mulitplying the expo- 
nential and additional constants. 



of a qubit, remains positive for all q > .5482 and N, 
which is equal to N3 since the results of the partial trace 
over any of the three qubits are equivalent, is non-zero 
for q > .2096. Once again the witness is successful in 
detecting tri-partite entanglement only up to a certain 
point despite the continued presence of tri-partite entan- 
glement. In addition, for .2096 < q < .5482 there is sur- 
viving tri-partite VF-type entanglement, as measured by 
N3, even when there is no longer any residual bi-partite 
entanglement. 

The effect of dephasing on the state p is to degrade the 
off-diagonal terms by e~ Kt . The state is detected by Wwi 
only for at < In 16— ln(— 5 + ^) while C, after tracing out 

any qubit, is max [| ( 4e~ Kt q - y/-3(q - l)(q + 3)),0]. 
We compare these values to N where again there is only 
one negative eigenvalue upon partial transpose of the 
density matrix: N = - min [^(3 - q(3 + 8\/2V Kt )), 0] . 
For q = 1 and sufficiently large Kt we do not find ESD 
for C or N but nevertheless find that the tri-partite en- 
tanglement cannot be detected by Wwi- For q < 1 ESD 
occurs for C and N as shown in Fig. fj] with the tri- 



IV. GHZ-TYPE STATE WITH BI-PARTITE 
ENTANGLEMENT 

We would like to explore a system that contains 
both bi-partite and tri-partite GHZ-type entanglement 
to compare and contrast the two types of entanglement 
evolution. We look at the following initial state 

1 — Q 

p=—^R + qpgb, (5) 

where p gb = \i)) gb {^\ and \^) gb = §(|000) + |111)) + 
i 1 1 10) . In order to detect the tri-partite entanglement 
of such states we must first find the proper GHZ and W- 
type entanglement witnesses. The closest state to \ip) gb 
with no W-type vectors is likely to be the GHZ state 
with an x- rotation about the third qubit. With this in 
mind we construct GHZ and W state witnesses Wc g b = 
§1 - p G (0) and W Wgb = ±1 - p G (6), for third qubit 
rotation of state pa by angle 8. The 9 that minimizes 
the witnesses is 8 ~ 28.075°. 
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With these tools we can explore the entanglement of 
this system. With no dephasing the state p is detected 
by yVogb for q > .763 and detected by Wwgb for q > -458. 
Taking a partial trace over one of the qubits leads to non- 
zero concurrence only when tracing over the third qubit. 
This concurrence, C12, is non-zero only when q > .529. 
The lowest lying eigenvalue of the partial transpose of 
p g b is when the partial transpose is taken with respect 
to the first or second qubit and gives non-zero negativity 
for q > .201. The tri-partite negativity in this case is not 
equal to the negativity and is nonzero only for q > .220. 

The effect of dephasing on the entanglement of this 
state is shown in Fig. [3] For q = 1 the tri-partite entan- 
glement witnesses cannot detect the entanglement below 
a certain threshold Kt, but C12, N,a,nd A3 do not ex- 
hibit ESD. For q < 1 all of the entanglement measures 
exhibit ESD. It is interesting to compare and contrast 
the behavior of the various entanglement measures espe- 
cially with respect to C12 ■ When q < 1 the entanglement 
that takes the longest amount of time to exhibit ESD is 
N. For sufficiently high values of q the tri-partite neg- 
ativity goes to zero before C\i. For the time between 
the subsequent sudden death of C12 and the ESD of N 
the entanglement present in the system is not of the bi- 
partite type since it is not destroyed by the partial trace, 
nor is it tri-partite in the sense that it can be measured 
by A3. The ability of the tri-partite entanglement wit- 
nesses to detect entanglement is weaker than any of the 
entanglement measures. 

As q decreases C12 decreases at a faster rate than any 
of the other entanglement measures. For sufficiently low 
values of q the concurrence, C12, exhibits ESD before A3 
such that the bi-partitc entanglement is the first to ex- 
perience ESD, followed by the tri-partite entanglement 
measured by A3. After this time there is entanglement 
that can be measured only by A. The behavior of A3 is 
different than that of the other entanglement measures. 
Rather than an exponential decay minus some constant 
there is initial exponential decay and then an inflection 
point before A3 goes to zero. Decreasing q futhcr we 
see that C\i undergoes ESD even before the W-state en- 
tanglement witness can no longer detect entanglement. 
This is in sharp contrast to the W-state case in which 
Tr(Wwip) always went to zero before the bi-partite C. 
This again shows the possibility of states with tri-partite 
W-typc cntanglmcnt without the presence of C when 
tracing over one of the qubits. 

The results of the above explorations are summarized 
111 the Table Q] There are a few points worth noting. 
First, we see the shortcomings of the entanglement wit- 
nesses in detecting mixed state entanglement. This is 
especially important for experimental realizations where 
entanglement witnesses may suggest there is no entan- 
glement present between systems when in fact there is. 
Second, for the GHZ and W-type states the different en- 
tanglement measures decrease at the same rate, i. e. , 
the exponential term is the same. However, this is not so 
for the state that has both GHZ and bi-partitc entanglc- 
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FIG. 3: (Color online) Demonstration of entanglement sudden 
death in p 9 (,-type states. Each plot shows C12 (solid line), N 
(large dashed line), A3 (medium dashed line), — Tr(WGgbPgb) 
(small dashed line), and — Tr(VVwgbPgb) (dotted line) as a 
function of the dephasing constant Kt. Top left: q = 1, the 
initial state is pure. In this case the tri-partite entanglement 
witnesses go to zero in finite time but the concurrence, neg- 
ativity, and tri-partite negativity do not undergo ESD. Top 
right: q = .9, for any sufficiently high value of q where q < 1 
the concurrence undergoes ESD before the negativity but af- 
ter the tri-partite negativity. Bottom: q — .6, the concurrence 
undergoes ESD before the tri-partite negativity and even be- 
fore the tri-partite entanglement witness. Thus, there is de- 
tectable tri-partite entanglement in the system though there 
is no concurrence remaining when one of the qubits is traced 



ment. Finally, we have found states where the negativity 
is non-zero despite the states having no tri-partite nega- 
tivity and no concurrence between any set of two qubits. 



APPLICATION TO QUANTUM ERROR 
CORRECTION 



With what we have learned concerning the entangle- 
ment behavior of tri-partite systems and the advent of 
ESD, we address the question of whether ESD affects the 
reliability of quantum error correction (QEC) or can ESD 
be used as a signature that error correction has only been 
successful up to a certain accuracy. Initial explorations 
of a possible connection were reported in Ref. [3 . Due 
to the fragility of certain quantum states QEC will be a 
necessary component of any working quantum computer. 
The multi-qubit logical states in which quantum informa- 
tion is encoded by a QEC code are typically highly en- 
tangled at the level of the physical qubits and thus may 
be subject to ESD. I explore the reliability of a qubit of 
quantum information encoded into a QEC code capable 
of correcting phase flips on one (physical) qubit. The 
goal is to determine whether ESD is a reliable signature 
to the success or failure of the code. 

To encode the state \tp) = a|0)+/3|l) into the phase flip 
code one simply uses the unencoded qubit as the control 
of two controlled-NOT gates each gate implemented be- 
tween itself and an additional qubit in state |0). This is 
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TABLE I: Summary of tri-partite entanglement decay behavior for GHZ-type, W-type, and states with GHZ and bipartite 
entanglement (GB) for expectation values of entanglement witnesses detecting tri-partite entanglement (the expectation value 
is negative when entanglement is detected), the negativity, the tri-partite negativity and the concurrence between the first and 
second qubits (all of which are the maximum between the values in the table or zero). 
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followed by Hadamard gates on all qubits. For simplicity 
we have chosen a,/3G$t. Unless a = 1 or (3 = 1 there is 
always some tri-partite entanglement present in this state 
which can be detected by a proper entanglement witness 
and measured by N (which in this case is equivalent to 
Ns)- We use the follwing W-state witness to detect all 
types of tri-partite entangement Wwh = |l — HpcH 
where H is a Hadamard on all three qubits. Detection of 
a possible error is done via syndrome measurement. 

When only one of the three physical qubits undergoes 
dcphasing there are parameters of the intial state and 
the dephasing parameter, p, which lead to finite time non- 
detection of entanglement in the error state p e \ by the 
entanglement witness Wwh- However, there is no ESD 
of N or N 3 (though N ^ N 3 ). The three qubit QEC code 
always corrects the error. Thus, there is no clear corre- 
lation between the detection of tri-partite entanglement 
by the entanglement witness or any of the entanglement 
measures and the success of the QEC. 

When all three qubits are placed in a dcphasing envi- 
ronment (where, for simplicity we have assumed that the 
dephasing parameter p is the same for all three qubits and 
thus ensure that N = N 3 ) the density matrix after the de- 
phasing, p e , exhibits tri-partite ESD, as shown in Fig. 2] 
Moreover, the QEC code will not, in general, completely 
protect the encoded quantum information. We can mea- 
sure the effectiveness of the QEC code by looking at the 
purity, P = Tr(pj) of the final single qubit state, pf after 
error correction and decoding or by looking at the fidelity 
of pf to the initial state. For the case at hand these two 
measures exhibit almost the exact same behavior. We 
can then ask whether the parameters for which this indi- 
cator is minimized (maximized) are those which exhibit 
(do not exhibit) ESD. We note that in general p/, and 
thus the purity, P, will depend on whether or not the 
syndrome measurement detects an error or not. The pu- 
rity of the final decoded state for both cases is shown in 
Fig.ffl 

When a = /3 = l/\/2 the QEC code perfectly cor- 
rects even complete dephasing on all three qubits depite 
the onset of tri-partite ESD. As a moves away from this 



point the purity decreases as a function of p, due to the 
inability of the QEC to correct the dephasing. Not sur- 
prisingly, some regions of low purity correspond to re- 
gions that exhibit ESD. However, even error regions in 
which the entanglement is detected by the entanglement 
witness can exhibit purity below ps .87. This suggests 
that, for the particular QEC code studied here, any cor- 
relation between ESD and the success of the code is not 
fundamental. A comparison of the onset of ESD due to 
the dephasing error with the purity is shown in Fig. [4] 



VI. CONCLUSIONS 

In conclusion, I have studied the effect of dcphasing 
on entanglement in three qubit systems. I have demon- 
strated ESD for both bi-partite and tri-partite entan- 
glement using concurrence, the negativity, and the tri- 
partite negativity. These were compared to the detection 
ability of tri-partite entanglement witnesses. In general 
we find that for mixed states the entanglement witnesses 
fail to detect the presence of entanglement well before 
the entanglement goes to zero. In addition, I have found 
a state in which there is no concurrence nor tri-partite 
negativity but entanglement as measured by the negativ- 
ity still exists. Based on these exploration I considered 
whether ESD affects the workings of the three-qubit QEC 
phase flip code and concluded that there is no fundamen- 
tal relationship between ESD and the failure of the code. 

As the number of system qubits grow studies of entan- 
glement become more complex as the number of different 
types of entanglement continually increase. However, a 
study of four qubit systems would allow for explorations 
of cluster state ent ang lement by using the proper entan- 
glement witnesses This would be of central impor- 
tance to issues of cluster state quantum computation [l9( . 

It is a pleasure to thank G. Gilbert and L. Vi- 
ola for helpful feedback and acknowledge support from 
the MITRE Technology Program under MTP grant 
#07MSR205. 
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FIG. 4: (Color online) Top left: Expectation value of entanglement witness — Tr(WwHPei) for single qubit dephasing error 
strength parameter p, before syndrome measurement and recovery, as a function of p and the initial state of the unencoded 
qubit parameterized by a. Lack of detection by the entanglement witness is evident despite the fact that error correction 
always works. Top center: — Tr(WwHPc) with dephasing of strength p on all three qubits. Lack of detection occurs at much 
lower values of p than when dephasing affects only one of the qubits. Top right: Negativity (which is equal to N3) of p e with 
dephasing of all qubits. ESD of negativity requires much higher values of p than lack of detection by the tri-partite entanglement 
witness. Bottom left: Purity, P, of final single qubit state after error correction and decoding when the syndrome measurement 
does not detect an error. Bottom center: P of final single qubit state after error correction and decoding when the syndrome 
measurement does detect an error. Bottom right: Comparison of P and the expectation value of the entanglement witnesses 
as a function of the initial state, paramaterized by a, and the dephasing strength, p. The contours are for — Tr(WwHPe) = -001 
(solid line), N — .001 (long dashed line), P = .99 if an error was not detected (medium dashed line), and P — .99 if an 
error was detected (small dashed line). The shape of the curves for the entanglement witness and N are similar though the 
entanglement witness fails to detect the entanglement significantly before ESD of N. However, while areas of higher purity 
and more entanglement are at lower values of p and values of a closer to l/y/2, there is no apparent correlation between the 
presence of entanglement and the success of the QEC code. 
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